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' Abstract 

, It is observed that interstitial hydrogen nucleii on a metalhc lattice are strongly 

|Zj I coupled to their near neighbours by the unscreened electromagnetic field mediating 

transitions between low- lying states. It is shown that the dominant interaction is of 
dipole-dipole character. By means of numerical calculations based upon published 
data, it is then shown that in stoichiometric PdD, in which essentially all interstitial 
^> ' sites are occupied by a deuteron, certain specific superpositions of many-site product 

. states exist that are lower in energy than the single-site ground state, suggesting the 

I existence of a new low temperature phase. Finally, the modified behaviour of the 

two-particle wavefunction at small separations is investigated and prelimary results 
suggesting a radical narrowing of the effective Coulomb barrier are presented. 
^ ■ 

. Key words: RDDI, phase transition, protons, deuterons, metal, interference, 

entanglement, fusion 
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1 Introduction 



Several metallic elements, notably palladium, vanadium, niobium and nickel, 
can reversibly absorb hydrogen up to the point of stoichiometry, in which every 
available interstitial site - of octahedral (O) or tetrahedral (T) symmetry - is 
occupied by a hydrogen nucleus. A single hydrogen nucleus in such an environ- 
ment exhibits a spectrum of singlet, doublet and triplet state representations 
of the local point symmetry group. The ground state is invariably a singlet 
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with even [+ + +] parity along each of the symmetry axes. The next level is, 

in an fee lattice, a triplet of states with parities [ — h +], [H h], [+ H — ] and 

an excitation energy of the order of 60 meV. The dipole moment between the 
ground state singlet and the first excited triplet is typically of the order of 
0.2Ae. Since the electronic Fermi gas couples weakly to the electromagnetic 
field quanta in this part of the spectrum, such a dipole moment gives rise to 
an essentially unscreened resonant dipole-dipole interaction (RDDI) between 
nearest neighbours [Ij. Simple geometrical considerations reveal this to be of 
the order of 20 meV per pair. Not only is this typically several times larger 
than screened static Coulomb interaction, it is also manifestly an appreciable 
fraction of the (on-diagonal) excitation energy of the dipole itself. Since in the 
quasi-stoichiometric loading regime each hydrogen has several nearest neigh- 
bours, it has previously been speculated [2] that there exist many-site states 
for which the total collective effect of the interaction is a multiple of the pair 
interaction. This paper sets out to answer the, in our opinion, intriguing ques- 
tion as to whether there exists any such collective state of quantum-entangled 
dipoles whose total energy is lower than the simple product of ground states, 
and to obtain an upper bound on the lowest possible energy of such an en- 
semble of coupled oscillators. 



2 Model 



The single particle states ipn are the solutions of 



2M. 



H 



■V^ + V{r) 



i!n{r) = e„^„(r) 



-where V(r) is the periodic potential experienced by an infinitely heavy posi- 
tive charge with fixed metal core positions. The reader is directed to [3] for a 
detailed discussion of the derivation of this potential using the DFT procedure. 
In view of the relatively large mass Mh of hydrogen nucleii, the lowest en- 
ergy solutions will generally be well-localised about local minima in V. These 
minima will coincide with sites of octahedral or tetrahedral symmetry in cu- 
bic lattices and hence the levels e„ are an assortment of singlets, doublets 
and triplet representations of the cubic point symmetry groups. The static 
(zero frequency) components of the potential disturbance due to the hydrogen 
nucleus is subject to a screening law of the approximate (Thomas-Fermi) form 



VHHir) 



-Kr 
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where K is proportional to the DOS at the Fermi surface. 
K is typically much greater than a reciprocal lattice vector. The static Coulomb 
H-H interaction between nearest neighbours is consequently small, typically 
not more than a few meV, and essentially state- independent. By contrast, the 
attenuation of the electromagnetic field due to a transition between levels is 
negligible over the dimensions of a lattice cell. 

Since the interparticle interaction is so strongly frequency dependent, the full 
Hamiltonian cannot be written in closed analytical form. However, matrix 
elements between pairs of two-site states are simply given by: 

g- A'|ri-r2-R|<5(ei J -ej^ )<5(ei2 -€j2 ) 

Hh,jv,i2,j2 = e^(ii,i2| , 5^ \ii,i2) (3) 

|ri — r2 — rl| 



where R is the intersite displacement vector. The 5(ej — e/) factors express the 
fact that only transitions between different levels give rise to an unscreened 
interaction. 

If just the lowest order term in the multipole expansion of the Coulomb op- 
erator I — 57 is retained, there is no need for double integration over both 

|ri — r2— R| ' ° 

volumes. In this approximation, ([3]) reduces to the familiar expression for a 
dipole-dipole interaction: 

Hh,ji;i2,j2 ~ ^(iikiNi) ■ 02^2^2) - ^(ji|R- ri|zi)02|R- r2|22) (4) 



If a classical dipole is located at every interstitial 0-site in an fee lattice, the 
interaction energy is lowest with the following orientations over a constant-z 
plaquette: 

< — o — > o < — o — > o < — 
oToioToio 

— > O < — O — > O < — O — !> (5) 

oioToioTo 
< — o — > o < — o — > o < — 



- where the open circles represent the sites of the metal cores at locations 
[100], [300], [111] etc. 

It can be shown that there is zero net interaction between parallel layers when 
each layer has such an arrangement. Guided by this classical analogue, we will 
limit our search for minimum energy states to those constructed from: 

[+ + +] parity states at all 0-sites in a 2; = plaquette 
[ — h +] parity states at 0-sites of even y and 
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[H h] states at 0-sites of odd y. 

[+ + +] ground states at all sites external to the plaquette 

For the rest of this paper we will use the shorthand \s^n) to denote the nth 
state of [+ + +] parity, \px, n) to denote the nth state of [ — h +] parity and 

\py, n) to denote the nth state of [H h] parity. The singlet ground state is 

accordingly written as |s, 0). If the site location needs to be made exphcit, we 
will append this in bold type thus: |s,0, 110). For clarification, we reproduce 
below an example of a pair of five-O-site states that are degenerate in zeroeth 
order and that are linked by the dipole-dipole interaction of (jl]): 



|s,0) o |s,0) 

O 1^3^,0) O 

|s,0) o |s,0) 



|s,0) o |s,0) 
o |s,0) o 

b.,0) o |s,o) 



or equivalently: 



|s, 0, 000) ® |s, 0, 020) O \py, 0, 110) ® |s, 0, 200) ® |s, 0, 220), 
0, 000) ® |s, 0, 020) (g) |s, 0, 110) ® |s, 0, 200) » |s, 0, 220) 



(6) 



The Hamiltonian matrix in the subspace of the two-O-site states 

|s, 0,000) ® 0,110) and 

\p^, 0, 000) (g) |s, 0, 110) is, according to (H]): 



H 



Cs.O + ^pfi 



^5,0 + ep,0 



(7) 



where the dipole length d = {s, 0\x\px, 0) = (s, 0\y\py, 0) 

The energy eigenvalues are in this case simply e^^o + ep,o ± ^^J''^^ 
More generally, and assuming full hydrogen occupancy and perfect lattice 
symmetry, the full many-site Hamiltonian matrix depends numerically upon 
just the lattice parameter a, the dipole lengths d and the energies e of the 
single-site states. 



3 Application to PdD 



Quasi-stoichiometric PdH and PdD are natural candidates for our model be- 
cause the adiabatic effective potential experienced by the hydrogen nucleus 
has been determined from ab initio DFT calculations [3M] and the theoretical 



4 



spectra found to agree well with IR spectroscopic measurements over a wide 
range of substoichiometric loading ratios. For this work we checked the results 
published in [3J by solving ([T]) on a real-space wedge mesh of pitch 0.03 A, us- 
ing the published effective adiabatic Pd-H potential and boundary conditions 
appropriate to the desired parity. The lattice parameter, corresponding to a 
displacement like [200] in our notation, is 4.07 A. The lowest few eigenvectors 
of the very large sparse matrix were solved using the dominant-diagonal iter- 
ation method, as in [5]. 

The lowest single-site energy levels, relative to the 0-site potential, were found 
to be: 



level 


PdH (meV) 


PdD (meV) 




82 


52 




151 


95 




233 


149 




289 


186 



- with the following dipole lengths along each of the three cartesian axes: 



dipole 


PdH (A) 


PdD (A) 


is 


0\x\p:,;, 


0) 


0.172 


0.153 


{s 




0) 


0.128 


0.112 


{s 


0\x\px, 


1) 


-0.003 


-0.002 


{s 


1 kba;, 


1) 


0.194 


0.165 



Double integration over both site volumes according to ([3]) yielded the follow- 
ing off-diagonal elements: 
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PdH (meV) 


PdD (meV) 


(s 


0,000 







llOlHlpy, 

\ \Ir y 7 


0,000) 




5,0,110) 


-27 


-21 


{s 


0,000 


(8) 





110|i/|p„, 


0,000) 




s, 1,110) 

7 7 / 


-20 


-15 


{s 


1,000 







llOliflpj., 


0, 000) 




5,1,110) 

7 7 / 


-15 


-11 


{s 


0,000 


® {Px 


1 


llOliflpj., 

1 l-t y 7 


0, 000) 




5,1,110) 

7 7 / 


-30 


-23 


{s 


0,000 


® (Pr 





200|i/|p^ 

1 Ix^-t 


0, 000) 


® 


5,0,200) 


-13 


-10 


{s 


0,000 


(g) 





200|//|p^ 

1 1 J 


0, 000) 


® 


5. 1,200) 


-9 


-7 


{s, 


1,000 


® {Px 


0, 


200|ii"b^, 


0, 000) 


® 


5,1,200) 

7 7 / 


-7 


-5 


(s, 


0,000 


® (Px 


1, 


200|i/|p^, 


0, 000) 


® 


5,1,200) 

7 7 / 


-14 


-11 


{s, 


0,000 


® {Px 


0, 


Q2Q\H\p^, 


0, 000) 


® 


5,0,020) 

7 7 / 


6 


5 


{s 


0,000 


(g) {pj. 





Q2Q\H\p.^ 


0, 000) 


® 


5.1,020) 


5 


4 


{s 


1,000 


® {Px 





020|i/|p^ 

1 1 J -t- 


0, 000) 


® 


5,1,020) 


4 


3 


{s, 


0,000 


(X) (Pr 


1, 


020|ii'|»^, 


0, 000) 


® 


5,1,020) 

7 7 / 


7 


5 


{s, 


0,000 


<8) {Px 


0, 


220|i/b„ 


0, 000) 


® 


5,0,220) 


-1 


-1 


{s, 


0,000 


(8) {p^ 


0, 


220|i/|;>„ 


0, 000) 


® 


5,1,220) 


-1 


-1 


{s, 


1,000 


® (px 





220\H\p^ 


0, 000) 


® 


5,1,220) 


-1 


-1 


{s, 


0,000 


® (Px 


1 


22Q\H\p,, 


0, 000) 


® 


5,1,220) 





-1 



3.1 5-site states 



As an illustrative example, we will calculate the lowest 5-site energy achievable 
for a given number of \p, 0) states. There are just five 5-site states having 1 
\p,0) and 4 |s,0) states, namely: 



|5, 0, 000) ® |s, 0, 020) ® |5, 0, 110) ® |5, 0, 200) ® \px, 0, 220) 
|5, 0, 000) ® |s, 0, 020) ® |5, 0, 110) ® \px, 0, 200) ® |s, 0, 220) 
|5, 0, 000) ® \s, 0, 020) \py, 0, 110) ® |5, 0, 200) ® |5, 0, 220) 
|5, 0, 000) (g) 0, 020) ® |5, 0, 110) (g) |5, 0, 200) g) |s, 0, 220) 
0, 000) |s, 0, 020) (g |5, 0, 110) |5, 0, 200) ® |s, 0, 220) 



(8) 
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The corresponding Hamiltonian matrix is (in meV): 



43 5 -21 -10 -1 
5 43 21 -1 -10 
H= -21 21 43 21 -21 
-10 -1 21 43 5 
-1 -10 -21 5 43 



(9) 



- where we have separated out the energy of the conventional ground state. 
The energy eigenvalues in this small subspace are 3,29,37,59 and 88 meV 
relative to Se^^o- For the 80 five-site states that comprise just one \p, 0) and 
four |s,n < 2) states, the lowest two energy eigenvalues are found to be -2 
and 26 meV relative to Se^^o- It is clear from this that the |s, 1) states make a 
significant contribution to the lowest energy eigenvector. 
In an attempt to find a converged value for the absolute mimimum site energy, 
Hamiltonian matrices were constructed for a series of plaquettes of increasing 
size up to a limit set by the memory capacity of our machine. The 15 sites 
included were: [000], [020], [110], [200], [220], [130], [310], [330], [420], 
[1-10], [-110], [240], [130], [040], [150] , in that order. 
The results obtained for the energies (i?o, Ei) of the lowest two states relative 
to iVe^^o are summarized in the following table. For the larger plaquettes, an 
energy cut-off was applied in order to limit the size of the matrix. 
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Sites 


p-states 


Cut-off (meV) 


States 


(meV) 


El (meV) 

J- \ / 


^0 / Site (meV) 


5 


1 




80 


-2 


26 





5 


2 




80 


31 


43 


6 


6 


1 




192 


-4 


19 


-1 


6 


2 




240 


12 


37 


2 


8 


1 




1024 


-15 


12 


-2 


8 


2 




1792 


-9 


12 


-1 


9 


1 




2304 


-16 


9 


-2 


9 


2 


300 


1044 


-14 


5 


-2 


12 


1 


300 


2784 


-21 


-7 


-2 


12 


2 


300 


3696 


-33 


-9 


-3 


12 


3 


300 


2200 


-34 


-12 


-3 


12 


4 


300 


4455 


-24 


-8 


-2 


15 


2 


300 


9660 


-40 


-23 


-3 


15 


3 


300 


5915 


-44 


-25 


-3 


15 


4 


250 


1365 


-26 


-8 


-2 



4 Conclusions 



The intrinsic complexity of this exact method and the inapphcabhty of a per- 
turbative approach have so far confounded our attempts to estabhsh a lower 
bound on the absolute minimum site energy. It follows from the variational 
principle that inclusion of higher |s, n) states, as well as further increase in pla- 
quette size, will result in even lower minimum energies. A mean-field approach 
is perhaps indicated, but we have as yet to find a sufficiently accurate formula- 
tion. It is nevertheless already clear from the above data that entangled states 
are favoured in the stoichiometric regime. The existence of a low temperature 
phase in which all the deuterons cohere in a mesoscopically entangled state is 
hence strongly indicated. 
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5 Over-cancellation of Coulomb barrier 



At small interparticle distances :- K\ri — r2\ << 1, the off-diagonal elements of 
the type we have been considering are comparable in magnitude, but opposite 
in sign, to the static Coulomb pair repulsion term. It is hence reasonable so 
suppose that, once coherence has been established, the height and width of 
the effective Coulomb barrier between neighbouring s,p state pairs is reduced, 
with a concommitant increase in the - normally infinitesimally slow - D-D 
fusion rate. In order to investigate the magnitude of this effect, we solved 
the two-particle Hamiltonian for the two states discussed in connection with 
d?]) above. Both the static ([2]) and dynamic interactions were included. 
Memory constraints limited us to a grid resolution of 0.06A. In view of ([6]) 
and ([7j), the solution was constrained to be of the form: 




where ijj is odd in its 5th argument and even in all others. 
It was found that the lowest energy solution, with 

e^e,,o + ep,o 5 (11) 



- exhibited an increased probability for close encounters of the two hydrogen 
nucleii right down to the limit of our resolution. At |ri — = 0.06A, the am- 
plitude was enhanced by about an order of magnitude over the simple product 
state that pertains when interaction is neglected. This exciting result implies 
that the dipole-dipole attraction effectively over-cancels the Coulomb repul- 
sion at least down to this length scale. The region of overlap was concentrated 
about the T-site lattice potential minima that are equidistant between the 
two 0-sites. 



6 Further work 



A search is currently being undertaken for other metallic lattices with high 
affinity for hydrogen and fiat effective potentials. A multi-level grid DFT al- 
gorithm of high accuracy has been developed for this purpose. 
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